Wavelet filters and infinite-dimensional unitary groups 

Ola Bratteli and Palle E. T. Jorgensen 



Abstract. In this paper, we study wavelet filters and their dependence on 
two numbers, the scale N and the genus g. We show that the wavelet filters, 
in the quadrature mirror case, have a harmonic analysis which is based on 
representations of the C*-algebra Om- A main tool in our analysis is the 
infinite-dimensional group of all maps T — > U (A'^) (where U (N) is the group 
of all unitary N-hj-N matrices), and we study the extension problem from 
low-pass filter to multiresolution filter using this group. 
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1. Introduction 

Digital filters are defined by a choice of coefficients which weight digital time 
signals from signal processing. It has been known since the eighties that these coef- 
ficients are also those which relate the operations of scaling and trans lation on th e 
real line M in the construction of multiresolution wavelets; see ( pill ) and ( ^l8| ). 
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There the digital filters are also known as quadrature mirror filters and they are fun- 
damental in generating wavelet bases in (R) . They can be determined by various 



rules which however at times appear somewhat ad hoc in the literature |JiMi91 



|Dau92[ , |BrJo97bt . Here we will describe a pairwise bijective correspondence 



between three sets: 

(i) the wavelet filters, 

(ii) a certain infinite-dimensional unitary group (also called a loop group), 
and finally 

(iii) a specific class of representations on Lp'iT) of some relations from C*-algebra 
theory. These last mentioned relations go under the name of the Cuntz 
relations and are widely studied in operator theory. 

The representations in ( ^ ) are not representations of groups, but rather repre- 
sentations of relations. These latter relations have been noted independently in the 
operator algebra community and in the subband filtering community; see especially 
Cun77| and [CMW92a|. In fact, the present authors stressed in |BrJo97b| and 



BrJo98[ | that the digital filters which go into the analysis and the synthesis of time 



signals with subbands are special cases of these representations. In the particular 
case of quadrature mirror filters (QMF) there are two bands, which are called high- 
pass/low-pass. In Section ||, we work out the general case of N bands, and make 
explicit the interplay between the separate viewpoints (tl)-(|iii|). This also has the 



advantage of explaining how wavelet packets (see |CMW92b]) arise very naturally 
from the representation-theoretic method in (0). Band filtering is described from 



an engineering viewpoint in |Vai93|. 

The question of reducibility versus irreducibility of the representations in (pl^) 
will be central, as will be unitary equivalence of pairs of representations. Our study 
of the equivalence question for representations is motivated by a desire for classifying 
wavelets and systematically constructing interesting examples. The question of 
reducibility of the representations in turn is motivated by the need for isolating 
when wavelet data is minimal in a suitable sense. Wavelets come with some specified 
numerical scaling which may be a natural number, TV = 2, 3, ... , or it may be an 
expansive integral matrix of some specified dimension d. (Here we shall restrict to 
d—\.) It turns out that the reducibility question for the representations is closely 
tied to whether or not the scaling data (in this case N) is minimal in a suitable 
sense; see Section ^ below. 

We also mention that the al gebra Om is used in the theory of superselection 
sectors in quantum field theory [DoRo89|. While Oat is used there in connection 
with representation theory, our viewpoint here is completely different: Doplichcr 
and Roberts use in [DoRo87| a category of endomorphisms oi On in deriving a 
noncommutative Pontryagin duality theory, also called Tannaka-Krcin theory, in 
giving an algebraic description of G when G is a given compact (n on-abclian ) Lie 

and 



DoRo89 



group. The intertwiners for systems of representations of G in 
|DoRo87 | induce the endomorphisms of O^. Our viewpoint is in a sense the 
opposite: We identify a class of representations of Oat which has the structure of 
a compact loop group, i.e., the group of all maps from T into U [N) where U [N) 
denotes the group of unitary N x N matrices. 
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2. Background on wavelet filters and extensions 

To fix notation, let us give a short rundown of tlie standard multiresolution 
wavelet analysis of scale N. We follow Section 10 in |BrJo97b], but see also 
|GrMa92| |, [ fVIeySTl ], | lMRV9e| ] and |Hor95| ]. Define scaling by N on L'^{R) by 

(2.1) {UO{x)^N-^^{N-'x), 
and translation by 1 on (R) by 

(2.2) (TO ix) - C(.T - 1). 

A scaling function is a function ip (R) , such that if Vq is the closed linear span 
of all translates T^(p, k G then ip has the following four properties: 

(2.3) {T'^ip ; A: G Z} is an orthonormal set in ^^(R); 

(2.4) Uip e Vo; 

(2.5) A„ezt^"Vo = {0}; 

(2.6) V„ezt^"^o-^'(M). 

The simplest example of a scaling function is the Haar function (p, which is the 
characte rist ic function of the interval [0, 1]. 
By (2.3) we may define an isometry 



(2.7) 

as follows: 

(2.8) 

Then 
(2.9) 

where ^ t-^ 
may define 

(2.10) 

and then 

(2.11) 

In particular, defining 
(2.12) 

we note that condition 

|2 



E„ bnfi ■ - n) 



m 



(2.13) PER(^|(^ 
which in turn implies 
(2.14) 



a-: 

- m{t)p{t) 

^ is the Fourier transform. If ^ G V-i = JJ^^Vo, then C/^ G Vo, so we 
= {UO e i'(T), 

mo{t) = m^{t) 
3|) is equivalent to 

k 

N-1 

\mo{t + 2TTk/N)f = N. 

k=0 



it) 



constancy a.e. t G 
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(Note that (2.14) does not conversely imply (2.13). Condition (2.14) merely im- 
plies that (p defines a tight frame. In the representation theory in Section |^ this 
distinction does not play a role. Only the unitarity of M{z) in (2.22) is relevant. 
See |H6r95, Theorem 3.3.6]. Similarly, the representation theory in Section Hap- 
plies in cases more general than the ones tied to the L^(R) wavelets via ( |j.l2[ ) 
and (2.18). In the latter case one also has the low-pass condition mo(l) ~ VN. 
In the loop representation ( 2.25 )-( 2.24 ), this is equivalent to Aoj(l) — 1/VN for 
j — 0, . . . ,N—1. These conditions involve the evaluation of the functions mo( • ) or 
Aoj{ • ) at z = 1. But z = e~^^ , where uj is the frequency variable, so the conditions 
are the low-pass conditions for w ^ 0.) 

If ^, 77 G U^^Vq, then we have the equivalence: 

£, _L T'^r] for all fc e Z 



Lfc=o "^c {t 

a- 



(2.15) 

If ^ e C/-IV0, then 
(2.16) 



t 



2TTk/N)mrj [t + 2T:k/N) = for almost aU t e 



— k), fc G Z, is an orthonormal set 
t 

Ef=o'K(i + 27rfc/A^)|'-iV. 

With the low-pass filter mg already given, now choose high-pass filters mi , . . . , rriN^i 
in L2(T) such that 

JV-l 

(2.17) {t + 2Trk/N)mj{t + 2nk/N) = Si^jN 

k=0 

for almost all i G M, and define wavelet functions ■01, • ■ ■ , '4'N-i by 

(2.18) VN^^{Nt) ^m,{t)ip{t). 
Using (|t|)-(|I|) one then shows that the functions 

(2.19) {T*-'^; ; fc G Z, i = 1,. .. , AT- 1} 

form an orthonormal basis for V-i n Vq ^ , and thus, using ( |2.5| )-( p^ , the set 

(2.20) {[/"T'^Vj ; n,fc G Z, i = 1,.. . ,iV- 1} 

forms an orthonormal basis for i^(M). Concretely, the functions in (2.20) are 

(2.21) U''T''i;^{x)^N~9ij,{N-''x-k). 



Thus, reformulating ( 2.17 ), orthonormality of {l/"T'°?/'i} is equivalent to unitarity 
of the matrix 



(2.22) M{z) 



moiz) 
mi (z) 



mo{pz) 
mi{pz) 



Too (p^ ^z'j 



\ 



mi 



N-l. 



for almost all z G T, where p = e 

N 



V TOAr_i(z) mN-i{pz) ■■■ mN-i{p^ ^z) J 

(so it is enough to consider z = e*^, < x < 



^). Here M is a function from T into U(Af) of a special kind. It will be convenient 
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for our analysis to consider general functions from T into U(A^). To this end, we 
do a Fourier analysis over the cyclic group Z/iVZ, and we introduce 

(2.23) ^,j(z) = i- V w-^m,{w), 



N 



and the inverse transform 

AT-l 

(2.24) m,(z) = -^'^^jO 

3=0 

These transforms are also the key to t he correspondences (|) <^ ([|) 
tioned in the Introduction. Then (2.24) may be summarized as 



4^ 



( 



(2.25) M{z)^A{z'"') 



1 

z 



1 

pz 



and the requirement that M (z) is unitary for all z G T is now equivalent to A(z) 
being unitary for all z e T. But A is an arbitrary loop. (A loop is by definition a 
f unct io n from T into U (A^). See the introduction to Section |[) Hence the formulas 
( ^.23| )-( 2.24 ) represent the bijection (§)-(|i|) alluded to in the Introduction above. 

One main problem considered in this paper is the problem whether the high- 
pass filters mi, . . . ,tojv-i selected such that (2.17) is valid can be chosen to be 
"nice" functions when the low-pass filter toq is "nice", where "nice" may mean 
continuous, differentiable, polynomial, etc. 

Using (2.23), this problem amounts to choosing a unitary matrix A{z) once the 
first row of A(z) is given as a unit row vector a{z). Ideally, one would like to use a 
selection map F from the unit sphere S*^^^^ of into A^-dimensional orthogonal 



P 



frames in 



(2.26) 



F(x) 



Fo(x) 



FAr_i(x) 



where the vectors Fo(x), . . . , Fjv-i(x) form an orthonormal basis for 
X G S*^^^^ and such that Fo(x) = x for each x, and define 

/ o.{z) 
Fi(a(z)) 

(2.27) A{z) 



for each 



V F^_i(a(z)) / 
For example, when N — 2 one may use Daubechies's choice 
(2.28) ¥i{xi,X2) = {~X2,xi). 

In general one can of course find measurable functions with these properties, but it 
is remarkable that one cannot choose them to be continuous when N > 2. Following 
Theorem 4.1 in [ [JiSh94|| and Re mark 10 .2 in [ |BrJo97bt , if S"""^ is the unit sphere 
in R", then a theorem of Adams [ Ada62| says that the highest number of pointwise 
linearly independent vector fields that can be defined on S"""^ is p{n) — 1, where 
the function p(n) is defined as follows: Let b be the multiplicity of 2 in the prime 
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decomposition of n; write 6 = c + 4c? where c G {0, 1, 2, 3} and d e {0, 1, 2, . . . }; 
and put p(n) = 2^= + 8d. Thus one verifies p{2) = 2, p(4) = 4, p(8) = 8, p(16) = 9, 
and in general p{n) < n \i n ^ {2,4,8}. Hence the only possibilities of finding a 
continuous selection function F are when = 2 or iV = 4. When N = 2 one 
can use Daubechies's selection function ( ^.28 ). When = 4 it is tempting to use 
the Cayley numbers to construct a selection function, but we will see that this is 
impossible. Note that if U(A'') is replaced by the real orthogonal group O(A^), a 
selection function can be found if and only if S {2,4,8} by using the complex 
numbers, the quaternions and the Cayley numbers, respectively. For example, for 
= 4, the matrix 



(2.29) 



is in 0(4) whenever the first row is a unit vector in 
C as 2 X 2 real matrices 





Zl 


Z2 


Z3 


Zi \ 




-Z2 


Zl 


~Z4 


zz 




-Z3 


Z4 


Zl 


-Z2 


V 


-Z4 


-Z3 


Z2 


Zl 1 



Using the identification of 



(2.30) 



Zl + iZ2 



Zl 
-Z2 



Z2 
Zl 



this corresponds to Daubechies's selection function ( 2.28 ). If A^ = 4, a se lection 
function does not exist in the complex case, by Proposition 3.6 in |]Jam76 ]. Spe- 
cializing to the case n = fc = A^ in this proposition, it says that, if the space of 
orthonormal A^- frames over admits a cross-section over the unit sphere in C^, 
and if p is any prime, then A^ is divisible by the smallest power of p exceeding 
(A^ - 1) / (p - 1). If Af = 4, this condition fails for p = 3 (while it is fulfilled for aU 
p if A^ = 2). (Note that a direct application of Adams's theorem only gives nonex- 
istence ofthecross-section for A^ > 4, as is also indicated in the proof of Theorem 
4.1 in |JiSh94].) (Note also that the discussion above settles the question implied 
by footnote 3 to the remark before Lemma 3.2.3 in iH6r95|l.) The claim made 



in | BrJo97b| , Remark 10.2] that there is such a selection function for A^ = 8 is 



definitely erroneous. What is true is that the matrix 



(2.31) 



v 

is in O (8) whenever the first row is a unit vector in R^. However, trying to convert 
this to an element of U (4) by the same trick as above leads to 



Zl 


Z2 


Z3 


Zi 


Zb 


Z6 


Zl 


Zfi 


\ 


Z2 


Zl 


-Zi 


Z3 


-za 


Z5 


zs 


-z-j 




z-i 


Zi 


Zl 


-Z2 


-Zl 


-zs 


Zb 


Zfi 




Zi 


-Zl, 


Z2 


Zl 


-Z8 


Zl 


-Zfi 


Zb 




Zb 


z& 


Zl 


z» 


Zl 


-Z2 


-Z3 


-Zi 




Zfi 


-Zh 


zb 


-Zl 


Z2 


Zl 


Zi 


-Z3 




Zl 


-Z8 


-z^ 


za 


Z3 


-Zi 


Zl 


Z2 




z% 


Z7 


-zq 


-zb 


Zi 


Z3 


-Z2 


Zl 


) 



/ 


Zl 


Z2 


Z3 


Zi \ 




-z^ 


'zi 


-Zi 


Z3 






Zi 


'zl 


-Z2 


V 


-Zi 


-'zi 


'zi 


Zl ) 



(2.32) 



but this is not in general unitary when the first row is a unit vector in C^. For an 
interesting short account of the matrices ( 2.29| )~( p.31 ), see [rau71|. The similar 
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selection problem in the context of splines and higer-dimensional spaces R" has been 



treated in RiSh91 and [RiSh92 ; see in particular pp. 142-145 in |RiSh91 



However, there is a completely different method of selecting A{z) using more 
global properties of the first row a{z) which applies in the case that a is a poly- 
nomial in z. This special case is interesting because it corresponds to compactly 
supported scaling functions ip: the scaling function tp constructed from mg by cas- 
cade approximation has support in [0, Ng— 1] if and only if the low-pass filter mo(z) 
is a polyn o mial of de g ree at mo st g — 1, and satisfi es too(1 ) = ^/N. See [ BrJo97b |, 



BrJo98f , |BrJo99|] , ||BEJOO|] , ||ReWe98| , and [t)au92|. In this context there is 



a method of extending a{z) to A{z) which is independent of the selection theory, 
and A{z) may even be taken to be a polynomial in z (with coefficients in B{C^)) of 
degree e q ual to t he degree of a. The embryonic form of this me thod is due t o Pollen 
Pol90|] , [lPol92[ , but has been developed further: see, e.g., ||ReWe98[ , |lVai93 |. 
We will give an account of the method in operator-theoretic form, and give the 



results in Proposition 3.2 and Theorem 4.1 
we therefore obtain 



Combining Theorem 4.1 with (2.25) 



Theorem 2.1. Let N e {2, 3, . . . } and let ip be a scaling function with support 
in [0, iVg — 1] where g G {1,2,3,...}. Then one can use multiresolution wavelet 
analysis to find associated wavelet functions ipi, . . . , V'w-i '^^^o having support in 
[0,7V.g-l]. 

Proof. The conditions on Lp, ■01, . . . , ipN-i correspond to the functions mp, 
mi, . . . , m7v_i being poly nomi als o f degree Ng 

and (|2.25|) 



follows from Theorem 
its proof. 



1, [H6r95 , hence Theorem 2.1 



More details are given in Corollary [4.2| and 

□ 



In Sections || and ^ we will consider the representations of the Cuntz algebra 
On defined by mQ,mi, . . . ,mjv_i when these are polynomials, and thus extend 
results from pEJOO[ . 



3. Decomposition of polynomial loops into linear factors 

We define a loop to be a continuous function from the circle T into the compact 
group U(iV) of unitary iV x iV matrices, where = 1, 2, . . . . The set C(T, U(A^)) of 
loops has a natural grou p structu re under pointwise multiplication, and this group 
is called the loop group [ PrSe86 |. We say that a loop is polynomial if its matrix 
elements are polynomials in the basic variable zgT = {wgC; \w\ = 1}. The set 
7^(T,U(iV)) of polynomial loops then forms a semigroup in C(T,\J{N)) called the 
polynomial loop semigroup. (Note that the formulas (|2.23D - (^.24|) for the bijection 
between the set of wavelet filters and the loop group are valid also in the wider 
category of L°°-functions , i.e., when each to^'^'' G L°°(T) and each Aij G L°°(T). 
This is the generality of | BrJo97b .) Note that if i7 is a given loo p, then U has 
a class in the Xi-group of C(T, Mn), i.e., Ki (C (T, Mn)) = Z, sec ||Bla86 



This 



class is called the McMillan degree in the wavelet literature | ReWe98| , and we will 



denote it by Ki (U). It can be computed as the winding number of the map 



(3.1) 



T 9 



det ([/(z)) G T. 



If U is in the loop semigroup, the integer Ki{U) is necessarily nonnegative. This 
follows from Lemma 3.1, below, apphed to u{z) — det {U{z)). The map U 
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Ki{U) is a group homomorphism from C(T, U(A^)) onto Z. We need an elementary 
lemma [ [ReWe98 |. 

Lemma 3.1. Ifw.T^Tisa polynomial, then u is a monomial. 

Proof. We have u = X]fe=o '^^z'^ for suitable coefficients Ck £ Z. By multiply- 
ing u by a nonpositive power of z we may assume cq 7^ 0. But then 1 = u{z)u{z) ~ 



CnCQz"+ powers of z between — (n — 1) and (n — 1) 
By induction, u{z) — cq. 



-CnCoz " and hence c„ 



0. 
□ 



We now use this to show that any polynomial loop U has a decomposition into 
linear factors | ReWe98 , p p. 60-6 1]. The idea in the proof of Lemma 3.3 goes at 
least back to | Pol90|] and | LLS96 , Lemma 3.3]. 

Proposition 3.2. Let U be a polynomial loop in V{T,\J{N)). The following 
conditions are equivalent. 



(3.2) 
(3.3) 



Ki{U) = d. 



There exist one- dimensional projections Pi 
V g U(^) such that 

U{z)^Ui{z)U2{z)---Ud{z)V, 

where 

U,(z) = zP, + {I - P,) 
for i — I, . . . , d and z G T. (// d = 0, then U{z) — V 



, Pd on C and a 



Proof. Since Ki{U^) = 1 and Ki{V) = 0, the imphcation (|3j) ^ (^ is 
obvious. We prove the other implication by using induction with respect to d. To 
this end, we will use the following reduction lemma. 



Lemma 3.3. Assume that 
(3.4) U{z)^Ao+Aiz + 



Akz' 



defines a U (z 7^(T, U(A^)), where A^ =/= and k > 1. Then U has a unique 
decomposition 

(3.5) U{z) = ((1 -Q) + zQ){Bo + Biz + ■ ■ ■ + Bk-iz^-^), 

where Q is the projection onto the range AkC^ of Ak {and then (1 — Q) + zQ and 
Bo + Biz + ■ ■ ■ + Bk-iz^^^ are m P(T, U(iV))). 



Proof. Let Q be the projection onto A^C . We have 



(3.6) 



Hence 
(3.7) 



1 = u{zyu{z) 

= {Al + A\z-^ + • ■ ■ + Alz-''){Ao +Aiz + - 
= AI.Aqz^'' + terms in higher powers of z. 



AlAo = 0, 



Akz'') 



and thus, since Q is the projection onto the range of A^, we obtain 
(3.8) QAa = and (1 - Q)Ak = 0. 
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Now put 

(3.9) W{z)^{{l-Q) + zQ)-^U{z) 

= {{1 - Q) + z-'Q)iAo + ■ ■ ■ + z'^Ak) 
= z'^QAo 

+ {il-Q)Ao + QAi) 

+ ■■■ 

+ z^{l — Q)Ak (ordered by increasing powers). 
It follows from (^^) that W{z) has the form 

(3.10) W{z) = Ba + Biz+--- + Bk-iz''-^, 



and ( |3.5|) follows. (In fact, = (1 - Q) Aq + QAi = + QAi, . .., and Bk-i = 
(1 — Q) Ak-i + Ak-) By our choice of Q, the decomposition ( |3.5[ ) is unique by 
(Fl). □ 



End of the proof of Proposition 3.2. We will prove the proposition by 
induction with respect to the McMillan index d. Assume first that d > I and that 
the proposition holds for all indices < cZ — 1 , and suppose that 

(3.11) U{z)^Ao + Aiz + ---+Akz'' 



has index d. But by Lemma 3.3, U has a decomposition 

(3.12) Uiz)^{{l-Q) + zQ)W{z), 
where Q is a nonzero projection. But then 

(3.13) d = Ki{U) = Xi((l - Q) + zQ) + Ki{W{z)) 

= dim{Q) + Ki{W{z)), 

so Ki{W{z)) < d — 1. (We define diniQ — rankQ when Q is a projection.) We 
may then apply the induction hypothesis to W(z). Finally, there exist dimQ one- 
dimensional projections Pi, . . . , -Pdim(Q) such that Q = Pi+P2 + - ■ ■ + Pdim{Q) , and 
then 

dim Q 

(3.14) ((1 -Q) + zQ)= II ((1 - P„) + zP„). 



It follows that U has the form (3.3). 

Finally, if Ki(U) — and U has the form ( ^.11 ), there are two possibilities: 



1. fc = and J7 is a constant unitary matrix. Then U already has the form in 



( |3.3D with d = 0. 

2. fc > 0. Then one may apply Lemma |3.3| to write t/(z) = {{I — Q) + zQ)W {z) 
with Q a nonzero projection, but as 

= KiiU) = dimQ + Ki{W) 

and Ki{W) > 0, this is impossible. 



This ends the proof of Proposition 3.2. The last argument in the induction chain 
is simpler if we use induction with respect to k rather than induction with respect 
to d. □ 
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4. Extensions of low-pass polynomial filters to high-pass filters 



In this section we will prove Theorem This theorem follows from Corollary 
4.2, below. Theorem 4.1 states that every vector a in which satisfies a certain 
orthogonality condition is the first row of coefficients of some element of the loop 
group. 

Theorem 4.1. Let a — {ao,ai, ... ,ag-i) be g row vectors in C'^ . The follow- 
ing conditions are equivalent. 

(4.1) The vectors satisfy the relations 

^ ("i I Oii+j) = for j ^ 0, 



("i l^i) = 1, 

i 

where we use the convention that a,; = if i < or i > g. 

(4.2) There exists a polynomial loop 

A{z) eV{T,\]{N)) 
of degree g — 1 such that the first row of A (z) is 

3-1 

i=0 

Proof. The implication ( |4.2| ) ^ ( |4.l| ) follows from considering the (0, 0) ma- 
trix entry of 

(4.3) A{z)A{z)*^l 

in Mn for all z, i.e., A (z) A (z)* is the constant Laurent polynomial 1. 

The other implication is proved by a very similar method as in the proof of 
Proposition 3.2. Again we use induction with respect to g. 

If g = 1, condition ( |4.l| ) just says that ag is a unit row vector, and we can find 
a constant function A{z) = A just by Gram-Schmidt orthogonalization. 

Assume that g > 1, and that the result has been proved for all smaller g. We 
may assume that Qfg_i ^ (otherwise we are already through by the induction 
hypothesis). Let P be the one-dimensional projection onto a^-i, i.e., 

(4.4) ag_iP = ag_i. 
Now define 

(4.5) P{z)^a{z){l-P + z-'^P) 

= {ao + aiz-\ h ag^iz^-^) {l - P + z^^P) 

= z^^aoP 

+ ao{l-P)+aiP 

+ z{---) 

+ ... 

+ zS-^ {ag^iP + ag.2 {1 - P)) 
+ zS-^ag^i{l-P). 
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But since (ao | ctg-i) = by (4.1), the z^^ term disappears, and (4.4) implies that 
the z^~^ term disappears. Therefore (3 {z) is a polynomial in z of degree g — 2. 
One now verifies from the unitarity of (l — P + z~^P) that the coefficient vectors 
of (3 satisfy the same relations (4.1) as those of a. Hence, applying the induction 
hypothesis, there exists a polynomial loop B (z) of degree g — 2 such that the first 
row of B (z) is (3 (z). Thus, putting 

(4.6) A{z) = B{z){l-P + zP), 



it follows from (4.5) that the first row of A (z) is a (z). This completes the induction. 



and thus the proof of Theorem 4.1 . 



□ 



Let us restate Theorem 2.1 and its proof in terms of the low-pass filter mo(z): 
Corollary 4.2. Let too be a polynomial, and let N £ {2, 3, . . . }. Suppose 



(4.7) 



|mo(z)| +\mo{zp)\ H h |too (zp^ ^) 



N 



for all z Cz T, where /? = e"" . Then there are polynomials {rrii ; i = 1, . . . , — 1} 
such that the combined system {mi ] < i < N} satisfies the unitarity property of 



In other words, every 



( p.22| ), or equivalently ( p. 17 ) with the convention z 
mo may be completed to a quadrature mirror filter system. 

Furthermore, when g £ N is chosen such that the degree of the polynomial toq 
is at most Ng — 1, then the polynomials mi, . . . , m^-i can be taken to have degree 
at most Ng — 1. 

Proof. With mo and N given as stated in the corollary, set 
(4.8) Ao^j (z) -^^ J2 ^"'"^0 ' ^ ^ ^■ 



N 



Then it foll ows from (4.7) that the coefficients aj, as a set of row vectors in 
satisfy (|]l|) in Theorem |4.il. Picking then A{z) = (A,j (z)) e 'P{T,V{N)) as in 
(4.2), using the theorem, and setting 



JV-l 



(4.9) 



m, (z) = ^'^^J- i^"") ' 



see ( [2.24 ) , we note that these functions satisfy the conclusion in the corollary. 

Finally, we note that, as in Remark 5.6, if too has degree at most Ng — 1, it 
follows from 



^) that Aqj (z) has degree at most g — 1] thus Aij (z) can be taken 
to have degree at most g — 1 by Th eorem 4.1, and then all to^ (z) have degree at 
most Ng — 1 by Remark 5_^ or (49). □ 

5. Representations associated with polynomial filters 

If Too, TOi, . . . , totv-i are complex functions on T such that the matrix A/ (z) in 
2^) is unitary, one may define isometries 5*0, . . . , Sn-i on (T) by 



(5.1) 
and then 
(5.2) 



— mi{w)^{w). 
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see [ |BrJo97b| , (1.16)-(1.17)] and, for iV = 2, |CMW92at (1.1)]. The isometries 
Sq, . . . , Sn-1 then have mutually orthogonal ranges, 



(5.3) 

and these ranges add up to 1, 
(5.4) 



s*s. = (5, J a, 



Af-l 



1. 



The relations ( p.3[ )-(5.4) are called the Cuntz relations (of order N), and the uni- 
versal C*-algebra generated by operators sqj • ■ ■ j sn-i satisfying these relations is 
called the Cuntz algebra (of order N). It is a simple antiliminal C*-algebra, which 
means that its space of Hilbert-space realizations cannot be equipped with a stan- 



dard Borel structure (see |Cun77|). In [BrJo97b| and Section |^ above, we identify 
intertwining operators W — JF* from L"' (T) onto resolution subspaces in (M) 
which intertwine the isometry So with the scaling operator / i-^ N^^^^f (x/N) on 
the appropriate resolution subspace in (M). The family of subspaces 



(5.5) 



■S^,L^ (T) 



/c = 0,1, 



e{0,i,...,iv-i}, 



corres pond s under W to wavelet packets (see [ CMW92b |) in (R). Relations 
(^^)-(5.4) make it immediately clear that the respective projections onto the sub- 
spaces (5^) are S.^^ ■ ■ ■ S^^ S*^ ■ ■ ■ 5* . In particular, these projections are clearly 



mutually orthogonal when k is fixed and multi-indices are varied 

Let us recall some results from 
space, and Si — > tt (s^) — Si €z B {TL) a representation of the Cuntz relations on 7i. 
Assume that there exists a finite-dimensional subspace JC dTC with the properties 



JKWOC] and [BEJOO . Let H be a Hilbert 



(5.6) 

and 

(5.7) 



S*IC c /C, 



0,1,. 



K is cyclic, i.e., [tt {On) IC] =H 



(here [7r(C'Ar)/C] denotes the closure of the linear span of all polynomials in 
So, ... , Sn-1 apphed to vectors in JC). Define operators Vi & B (JC) by 



(5.8) 
Then 

(5.9) 

Define a map a on B {IC) by 
(5.10) 



V* 



N-l 



E ^^^^ = 1^ 



i=0 



N-l 



a{A)=J2v,AV*. 



1=0 



Then cr is a unital completely positive map, and 



Theorem 5.1. ( pJKWOOt ) There is a positive norm-preserving linear iso- 
morphism between the commutant algebra 

(5.11) TT {On)' = {AeB{n); At: {x) ^ it {x) A for all x € On} 
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and the fixed-point set 

(5.12) B {ICf = {A e B {IC) ; a (A) = A} 
given by 

(5.13) TT (On)' 3 a — y PAP, 

where P is the projection oJTL onto IC. In particular, tt is irreducible if and only if 
a is ergodic. 

More generally, if ICi, IC2 {with corresponding projections P*-^-* and P*^^-*) are 
S* -invariant cyclic subspaces for two representations tti, 772 of On on Hi, Ti.2, and 

(5.14) T/«^pOV^.(.,)|^^ 
for j — 1,2, i = 0, . . . , N — 1, define p on B (/Ci, /C2) by 

(5.15) p(A) = ^v;(^)Ai^w*. 

i 

Then there is an isometric linear isomorphism between the set of intertwiners 

(5.16) {Ac,B{rLi,rL2) ; Atti (x) ^ 1^2 {x) A for all x e On} 
and the fixed-point set 

(5.17) {PeS(/Ci,/C2);p(B) = B} 
given by 

(5.18) A — >P = P(2)AP(1). 

We argued in BEJOOj , (4.14)-(4.18)] that if mo, . . . ,mN-i are polynomials 
in the circle variable z, N = 2, and H = (T), then such a finite-dimensional 
subspace IC exists, having dimension Ng and spanned by 1, z~^,z~^, . . . , z~^s^^. 
In this section we will extend this result to > 2, and we will see in Proposition 5.5 
that we can do slightly better than what the N = 2 result indicates. Let us assume 
from now on that the loop group element A{z) in (2.2S)-( 2.25| ) is a polynomial 
of degree g — 1, so that mo (z) , . . . ,r7ijv-i (z) are polynomials of degree at most 
N {g -1) + N ~1 ^ Ng -1. 

We will use the notation A{z) — {A^j {z))fjj,^ for the loop-group element 
A: T ^ \J (N). Since the Fourier expansion is finite, A (z) has the form 

(5.19) A(z) = ^z'=AW, 

k=Q 

where A^'''^ G B (C^) for /c = 0, . . . , g—l- The factorization in Lemma [s!^ motivates 
the name genus for g. 
We have 
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Corollary 5.2. // A (z) is a general polynomial of z with values in B (C"^) 
of the form ( [S.1S| ), the following four conditions (5.20)-( ^.23| ) are equivalent: 

(5.20) A (z)* A (z) = In, z eT, i.e., A takes values in U (A^); 



(5.21) = 
(5.22) 



and 
(5.23) 



with the convention that 



Ijv ifn = 0, 
ifneZ\{0}, 
A(")==0z/m^{0,l,...,5-l}; 

there are projections Pi, . . . ,Ps in B (C^), positive integers 
ri, . . . ,rs, and a unitary G U (N) such that 

A {z) = (n;=i (Ia' - P, + z'^Pj)) W; 



there are projections Qi,Q2, ■ ■ ■ , Qg-i o.'^d o- unitary V € V (N) 
such that 

9-1 

A^^^ =vJ2i^N-Ql)■■■ 
j=^ 

■ ■ ■ (iw — Qj-i) Qj (iw ^ Qj+i) • ■ ■ 

• • • (liV — Qg-l) , 



S-1 



Proof. This follows from Proposition B.2 and explicit computations 



□ 



Remark 5.3. The case g = 2 = N includes the family of wavelets introduced 



by Daubechics |Dau92| and studied further in |BEJO0]. Note that g ~ 2 yields 
the representation 



(5.24) 



by ( 5.25 ). But then ( 5.21 ) takes the form 
(5.25) 1a, - Q, 

which will be used in the sequel. 



^(i)*^(i) ^ 



Let us return to the connection between loop-group elements and representa- 
tions of Ojv- The algebra On has the following basic representation Si Si on 
(T): 

(5.26) Sj^ [z) ^ z^i (z^) , < j < TV, e e (T) , z e T. 



If Si Ti is any other representation, then as noted in |Cun77|, |Cun80 , and 
|BJP96 |, the system {S*T^ . . realizes a unitary N x N matrix with entries in 

6(l2(T)). 
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Proposition 5.4. Let (Si) be the basic representation of On, and let (Ti) be 
an arbitrary representation. Then the following two conditions are equivalent: 



(5.27) 
and 
(5.28) 



Each operator S*T^ on (T) is a multiplication operator; 



the Ti-representation has the form 

(z) = m, (z) e (z^) , 0<i<N, ^GL^T), zGT, 

where thq, . . . ,tojv-i € L°° (T). 

Proof. (|^ ^ (|]28|) : If (^f^) holds, then there are functions ^ e L°° (T) 
such that S*T^ = Aij where the right-hand side also denotes the multiplication op- 
erator determined by the function in question. Then 

N-l N-1 N-1 



(5.29) 

Setting 

(5.30) 



— ^ SjS*T,-^ - ^ SjA^j - ^ A.,^^ 

J=0 j=0 3=0 



N-1 

mi (z) = Ai^j 



and using the Cuntz relations, we conclude that Ti has the form in (5.28). 

(|5.28|) ^ ([3-27|) : If, conversely, the representation (Ti) is given to have the form 
in ( p. 28 ), then one checks that the filter functions must satisfy the unitarity 
property ( 2.22] ) above. If then Ai j (z) are given by then 



(5.31) 



and we conclude that (5.27) holds, i.e., S*T^ is multiplication by the function 
Aij ( • ), where Aij is derived from m.^ via (2.23) □ 

We will henceforth denote the representation Ti defined by ( p3.31 ) by Tj;'^\ so 
in particular. 



T, 



(1) 



(5.32) 

where the notation Si is reserved for t he r e presen tation defined by ( ^.26| ). 

We will now apply the results in (5.6)-( 5T^ ) to analyze these representations 
further. As mentioned after ( 5.18 ), the linear span K, oil , z~\ z~^, . . . , z~'^ is cyclic 
and r*-invar iant for a suitable r € N, i.e., satisfies (5.6)-(5.7) with T in lieu of S", 
see [ BEJOO , Proposition 3.1 and Corollary 3.3]. If m p, . . . ,mN-i can be derived 
from a polynomial loop A (z) of degree g — 1 by ( 2.24 ), we may explicitly estimate 
r. To this end, let us look at the action of T^^^ * on e„ (z) = z". Put n = j + Nl 
uniquely, where j S {0, 1, . . . , iV — 1} and I £ Z. Then, using (5.29), 

(5.33) i;^^^ *e„ = S*^A~^TP^)en 



EE^ 



-Nk- 
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ei-k if j' = j, 



But n- Nk ^ j + Nl - Nk ^ j + N {I - k), and since 

(5-34) S*,ej+Nii-k) = |q 

we obtain 
(5.35) 



ifjVj, 



9-1. 



7n('4)* _ .(fc) 



k=0 



(A) + 

It follows that Tj is represented by a matrix which is a slanted block matrix of 
the following form: 



(5.36) 




The matrix elements outside the shaded blocks are all zero, and each block has N 
columns and g rows. Each block is a translate N steps to the right and 1 step up 
compared to the previous one, and the most central block is located with corners at 
(0, 0), {N - 1, 0), {N -l,-g + 1), (0, -g + 1). It follows that if we fix one n e Z, 

(A) * 

iterated applications of for various i's ultimately will transform this vector 

into a linear combination of em's where m € {0, —1, —2, . . . , — r}, where r is the 
largest integer such that (— r, — r) is contained in one of the blocks. Thus the linear 
span /C of {eo, e_i, . . . , e_r} is Tj*-"^-* *-invariant, and because of the relation 



(5.37) 



1 



iA)rpiA)* 



IC will also be cyclic (see [ BEJOC , Section 3] for more details of this argument). 
Thus, computing r more explicitly, we deduce 

Proposition 5.5. IfT^^^ is the representation of the Cuntz algebra On de- 
fined by a polynomial loop A (z) of degree g — 1, then the subspace 

(5.38) /C — linspanjeo, e_i, e_2, • • ■ , e-r} 



satisfies (5_^) and (5/7), i.e., K. is T'^'^^ * -invariant and cyclic, where 
(5.39) r = g- 



5-1 




gN-1 


_N - 1_ 




_N -1 _ 
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Here \_x\ is the largest integer < x. 



Proof. From the figure (5.36) it follows that 



(5.40) r = 




if g < N, 

if AT + 1 < .9 + 

if 27V + 1 < 5 



1 < 2iV, 
h 2 < 3A^, 



I.e. 
i.e. 
i.e. 



g-KN - 
N ~\<g- 
2N -2<g- 



1, 

1<2N -2, 
- 1 < 37V - 3, 



Here the first column of ranges ("if . . . ") is taken from the box diagram ( 5.36 ), and 
the second column of ranges ("i.e., . . .") is derived from the first by subtraction of 

i- arises in (5.39). □ 



suitable integers. The second column shows how the ratio 



Remark 5.6. Proposition 5.5 could alternatively have been proved by exactly 



the method used to show (4.18) in |BEJO0] from Remark 3.2 there, using joint 
invariant sets for the maps n ^ {n~ k) /N for fc = 0, 1, . . . , gN — 1. Note that a 



consequence of (2.23) and (2.24) is that the degree of all the polynomials Ai,j (z) 
for j, j = 0, . . . , iV — 1 is at most g — 1 if and only if the degree of all the polynomials 
rrii (z) for i = 0, . . . , iV — 1 is at most gN — 1. 



The case considered in detail in [BEJOOj was TV = .g = 2, so we regain the 
result r = 3 from there, i.e ., in that case K, — span {cq, e_i, e_2, 6-3}. Note that 
the denominator — 1 in ( 5.39 ) shows the dependence of r = (dim/C) — 1 on the 
scaling n umbe r N for the given wavelet filter. So when N > 2, the d enomina tor 
TV — 1 in ( 5.39| ) yields a smaller value for r than the value gN — 1 from | BEJO0|| in 
the special case N = 2. 

We will now apply Theorem |5.l| to t he representations above, when K. = JC^"^^ 
is given by ( 5.38 ) and r — r'-^^ by ( |5.39 ). The operators Vi & B (/C*^'^-') defined by 

in lieu of Si in (5.8) will be denoted by V^^\ and then the a defined by (5.10) 



is denoted by a'^'^^ 
(5.41) 



Xa,a) 



N-1 

E 

j=0 



for X e 6(/C(^)). If r(^), r(^) are two representations of our kind, let a^^-^^ 
denote the p defined by ( 5.15| ), i.e., 

(5.42) a(^'^Hx) = Y.v} 



so o-(^'"^) is in B (/C'^'^^/C*^^'). We will mostly be interested in the situation that A 
and B have the same genus g, and then /C*^'^^ = /C*-^^ by Proposition 5.5 (of course 



we may always replace g'"^', g''^^ by max {5'^"*^ g'^^ } and thus assume /C'-^-' = 
/C(^)). Now B {K:'-^\ K.'-^^) can be made into a Hilbert space K.'-^-^^ in a natural 
fashion by defining the norm of V G S(/C(^),/C(^)) by 

(5.43) UWi^'^'' = Tr^-*^ (i/;» = Tr^-^) (iP^j;*) , 

where Tr*-'^-', Tr'^-* are the standard unnormalized traces on B^K,^'^^), B{lC^^^), 
respectively. Then a^^'^^ e B (/C'-^''^^). Because of Theorem our main con- 
cern will be the spectral multiplicity of 1 in the spectrum Sp (cr^^^"^)), as well as 
the associated eigensubspace IC^^'^K It turns out that these items are much sim- 
pler to compute for the adjoint cr^^"^)* G B (/C^-^'^^) (i.e., the adjoint of cr(^'^) 



18 



OLA BRATTELI AND PALLE E. T. JORGENSEN 



as a Hilbert-space operator). Then Sp (cr'^'"^'*) — Sp (cr^^'"*)), and in particular 
1 G Sp (cr(-^'^)) if and only if 1 G Sp (cr^-^^^)*). Vi A = B and the completely posi- 
tive unital map ct("^''^) admits a faithful invariant state, it follows from | BrJo97a , 
Lemma 6.3] that the dimensions of the fixed-point sets of cr^"*'"^^ and g'^^'^^* are 
the same. However, we will see that cr'-"*''^' does not necessarily admit a faithful 
invariant state. Our strategy will be to compute Sp (cr'-^''^') by evaluating 



(A.B): 



det xl-ai^'B)) 



(5.44) det 
and then compute the fixed points by hand. Note that from (5.42) we have 

(5.45) a^^'^^*{X) = Y,vl^^*Xv}^^ 

i 

for X e IC'^^'^^ = B {IC^^\ K.^^^) . If X = E_k,-i = the rank-one partial isometry 
mapping e_i into e_fc, we compute, in Dirac's notation, 

(5.46) (i?_,,_,) ^Y.^'''E->^:~i^^*- E k/''^--^ 



which is complicated to handle. On the other hand, 
(5.47) a(^^^)*(£;_fc,_0=E^/ 



-EN 



B) * 



e-k 



-EN 



B) * 



T- e-i 



where the last step used T*-invariance of /C, see ( |5.6D and (p. 81). But T* is given by 
( 5.35 ), and hence the matrix elements of (t(^^^) can be computed explicitly. Let us 
do the calculation in a special case: 



The case g = 2 and N > g. Then r = 2 by Proposition 5.5, and K, is 
three-dimensional. Then the loop A (z) has the form 

(5.48) A{z)^Vil-Q + zQ) = A^°^ + zA^^^ , 

where 1^ is a unitary and Q a projection in B (C^). As noted in ( 5.25| ), we have 

(5.49) A^°^*A^°'> = 1n^Q, A^^'>*A^^'> =Q. 
Now define 



(5.50) 



A,,; 



(1 - Qk, = s^.j 



Qi,j 



z,j = 0,l,2,...,iV-l. 



Setting Xi := Xi,i, then ( 5.47 ) and ( 5.35| ) give 

a(-4'^) * (Eo.o) = Ao-Bo.o + (1 - Aq) -B-i.-i, 

^{A,A) * = AAr_i_E_i^_i + (1 — Aat-i) £'_2,-2, 



(5.51) 



y{A.A)' 



{E-2,-2) — Xn-2E-1^-1 + (1 ^ AAr_2) E^2-2, 



CT^"'"'" (-Eo,-l) — Ao,jV-l£'o,-i — Ao,Ar-i£'-i,-2, 



JA.A)* 



{Eo^-2) — Ao,JV-2£'0,-l — Ao,JV-2£'_i^_2, 



CT^"'"'" (^^-1,-2) — AAr_i^Ar_2-E-l,-l — AjV- 1, ^-2-^-2,-2 , 



yiA.A)* 
M^^)*{E^k,- 



^(A,A) * {E^i^^kT whenever Q<l<k<2. 
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From here one computes 

(5.52) det (xl - a^^-'^A = dot (xl - a^^^^^ 



^ X {x-l){x- Ao) (a; - {Xn-i - {x - Xo,n-i) {x - Xo,n-i) ■ 

Since the Xij are matrix elements of a projection, we see immediately that the 
eigenvalue 1 has multiplicity one except in special cases as when Aq = 1, AAr_i = 1 
and AAr-2 = 0, or Ao,jv-i = 1- Hence: 

Theorem 5.7. If g — 2 and N > 2, the representation T^"^^ is irreducible for 
generic loops A. It is reducible if and only if the projection Q ^ B (C^) defining A 
by ( 5.48| ) satisfies at least one of the conditions 

(5.53) go,o = ( ^ Ao = 1) 



(5.54) Qn-im-i = and Qn-2,n-2 = 1 
Proof. We will use the fact that the A 



( <;=^ Aat-i = 1 and A7v-2 = 0). 
numbers are the matrix entries of a 



projection, 1 — Q. Since (1 — Q) = (1 — Q), we have 



(5.55) 



At — Ai, 



Thus the case ( 5.53 ) may occur with 



(5.56) 



Q = 



V 






••• \ 




















Qicd 












= (0) e grcd- 



Also ( p3.54 ) may occur with 



(5.57) 



Q = 












^ 
















Qrcd 























••• 





1 








••• 












Qrcd ffi 



1 




But the case Qq 



T, i.e., A 



0,Af-l 



A 



N-1,0 



Z, since Q is. Specifically, 



1, cannot occur. To see this, 
recall that, by (5.25) and (5.5C), the Xij numbers are the matrix entries of Hat — Q, 
which is a projection in 

(5.58) 

where {ejj^Q^ is the canonical basis for C''^. So we are considering the possibility 
of having an off-diagonal entry Ao,7v-i assume the value 1, and that is impossible. 
In fact, let P = P* = be any projection, and let e, e' be given orthogonal unit 
vectors. Then 



(5.59) 



MPs') 



1 

^2' 
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and so in particular the value 1 is excluded. Note that ( ^.59 ) is clearly sharp: take 
'1/2 1/2^ 



1/2 1/2 
(5.60) 



To prove ( p.59D , note that 

\{e\Pe')\' + \{e'\Pe')\^<\\Pe'\\\ 



by Bessel. Since {e'\Pe') = llPe'll , we get 



(5.61) |(e|P£')l < Vll^^'ll 11^^11= ll^e'll- l-ll^^' 



-'l|2 



1/2 



1 



\Pe'\\ • 11(1 - P) e'W < \ (\\Pe'\\' + ||(1 - P) e'\\') = - 



2 _ 1 

^ 2' 



This proves (|5.59|) 



□ 



Remark 5.8. So in summary, we m ay de fine a transformation a for any system 
of n umbe rs A^j- via ( 5.51 ), or the matrix ( |5.63| ) below, and its spectrum will be given 
by ( 5.52 ). We arc then interested in when points from the spectrum of a can attain 
the value 1, and we saw that some can, and others cannot, viz., if the A^j are 
derived from an A with g = 2, then we get the following estimates on the points 
Ao, Aat-i — AAr_2, Ao,iv-i, and Xn-i,o in the spectrum of a: 

0<Ao<l, -1 < Aat^i - AAr_2 < 1, 

and, by the middle step in ( 5.6l|) , 



|Ao,jv-i| = IAat-i^oI < min {(Ao • (1 - Aq))'/' , (A^^i • (1 - Xn^i)^^} < ^. 

In particular, the spectral radius of such a, a — cr''"*-' is at most 1. Nonetheless, 
we know that cr^'^-' is generally not contractive in the Hilbert space B (/C^'*'') of 
Hilbert-Schmidt operators, although cr^'^'' obviously is contractive as a completely 
positive map on the C*-algebra B (/C^"*^). 

Remark 5.9. Hence for a given unitary V GlJ {N) and projection Q E B (C^) , 
we see that the eigenvalue 1 may have the following multiplicities in the character- 
istic polynomial for various values oi A ^ {V, Q): 

1 has multiplicity 3: This occurs in one case. 



Case 1. Ao = 1, Xn-2 ~ 0, and Aat-i = 1. The projection 1 — Q then has the 



form 



(5.62) 



1 





... 

























P : 

























... 














... 





1 . 



where P is a projection in B (C^ ■^) (so when = 3 there is no more choice. 

1 has multiplicity 2: This occurs in two mutually exclusive cases. 
Case 2. Ao = 1 and {Xn-2 > or Aat^i < 1). 
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Case 3. Ao < 1, Xn-2 = 0, and Aat^i = 1. 
1 has multiplicity 1: This occurs in all remaining cases. 

Since the dimension of the fixed-point set of a^^^ is always at most equ al to the 
multiplicity of 1 in the characteristic polynomial, it follows from Remark p.9| and 
Theorem 5.1 that the linear dimension of the commutant of the T^^\= T^''^-')- 
representation is always at most 3. Since the smallest-dimensional nonabelian C*- 
algebra is B (C^), which has dimension 4, it follows that the commutant is always 
abelian when g — 2, and the representation decomposes into a sum of at most three 
irreducible mutually disjoint representations. In order to get the exact number, we 
have to compute the fixed-point space of fj*-"^' exactly. To this end, we note that the 
matrix of cr'^'^\ relative to the basis {Ei k ; —2 < i, k < 0} of B (C'^) is (we actually 
compute the matrix of a^^^ * using ( |5.5l| ) and take the adjoint) 
(5.63) 





(0,0) 


(-1,-1) 


(-2.-2) 


(0,- 


1) 


(0,-2) 


(-1,-2) 


(-2,-1) 


(-1,0) 


-2,0) 


(0,0) 




1-Ao 

























(-1,-1) 







1 — Aat-I 






















(-2,-2) 







1— A^_2 






















(0,-1) 













-1,0 





— Ajv-1.0 











(0,-2) 













-2,0 





— Ajv-2,n 











(-1,-2) 







— Ajv-2,W-1 






















(-2,-1) 







— AjV-l,W-2 






















(-1,0) 






















— Ao,Af-l 


Ao,Af-l 





(-2,0) 






















— Ao,Af-2 


Ao,Af-2 






Thus one may compute the eigenvectors of cr'-"^^ . The generic result is as in Table 
i 

Table 1. Eigenvalues, multiplicities, and eigenvectors of cr^'^^ 



Multiplicity 


Eigenvalue 


Eigenvectors 


4 





Eo,-2, E-2fl, -Eo.-l + -B-1,-2, E-i^o + E-2,-1 


1 


Ao 


Eo,o 


1 


1 


1/C = Eq^o + E-i,-i + E-2,-2 


1 


Aat-I — AAf_2 


(1- Ao)(l-AAr_i)£;o,0 

+ (AtV-I — A7V-2 — Aq) (1 — Aat-i) 
— AAr_2 (AjV-1 — AjV-2 — Ao) -E-2,-2 


1 


Aat-I^o 


Aat-I^o^-O.-I + AAr_2,o£'0,-2 


1 


Ao,Af-l = Aat-i^o 


Ao.Af-l-E-l.O + Ao,JV-2-E'-2,0 



Note that in non-generic cases, like Ao — 1, and/or (Ajv-i = 1 and Xn-2 — 0), 
some of the listed eigenvectors are zero. This is not surprising since these are 
exactly the cases where the root 1 in the characteristic polynomial is multiple. Let 
us consider these cases separately: 
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Case 2. Ao = 1 Ao,jv-i — Ao,jv-2 — 0) and Aat^i < 1 or \n-2 > 0: In this 
case the last two eigenvectors in the list in Table |] are zero, but one may verify 
that these may be replaced by the vectors -Eo,-2 and £'-2,0- Thus the eigenspace 
corresponding to 1 is indeed 2-dimensional in this case, and spanned by l^c and 

-£'0.0- 



Case 3. Ao < 1, Ajv_i = 1, Ajv-2 = 0: In this case it follows from ( 5.55 ) that 



^N-i.k = for aU fc 7^ — 1 and Ajv-2,fc = for all fc 7^ TV — 2, so the matrix in 



(5.63) degenerates into the upper-left-hand 3x3 matrix, which is 

/ Ao 1 - Ao 

(5.64) 10 

\ 1 

So we see that a new eigenvector for Aat-i — Ajv-2 = 1 is 

(5.65) E^2,-2. 

Let us now look at the even more degenerate case where the multiplicity of 1 
in the characteristic polynomial is 3: 

Case 1. Ao = 1, AAr_i = 1, AAr_2 = 0: From the discussion of Case ^ above 



it follows that then the upper-left 3x3 matrix of (5.63) is I3, while the rest of the 
matrix is zero. Thus the fixed-point set is the linear span of 

(5.66) Eo^o, and E_2-2- 

Thus the dimension of the eigensubspace is in all cases where the parameters 



are restricted as in Remark 5.8 equal to the multiplicity of the root 1 in the char- 
acteristic polynomial, and the above analysis implies 

Corollary 5.10. If g — 2 and N > 2, the representation T*-^^ decomposes 
into at most three irreducible representations which are mutually nonequivalent. 
Generically, T'^^^ is irreducible, and otherwise the decomposition structure is sum- 
marized in Table ^ 



Table 2. Decompositions from eigenvalue-one multiplicities, and 
eigenvectors for a^^^ 



Case 


Number of 


The fixed-point set of cr^^) 




subrepresen- 


is spanned by 




tations 




Ao = 1, AAr-2 = 0, 


3 


£'0,0, £-1,-1, and E-2,-2 


and Aat-i ~ 1 






Ao < 1, Aat-i 1, 


2 


Ik and E-2.-2 


and AAr-2 = 




(i.e., £0,0 + £-1.-1 and £-2,-2) 


Ao — 1 and 


2 


Ik and Eq^q 


(Aat-I < 1 or AAr_2 > 0) 




(i.e., £0,0 and £_i _i + £-2-2) 



Remark 5.11. Table |^ shows a remarkable feature which these representations 
share with the g = 2, N = 2 representations studied in |BEJOO|, namely that the 
fixed-point set of a(^) is both abelian and an algebra in all cases. Hence it follows 
from (5.13) that the projection P onto /C commutes with the commutant n {On)' , 
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and thus we may find cyclic vectors for the various subrepresentations by picking 

vectors in the range of the minimal projections in the fixed-point set B (/C) . For 
example, for the first case in Table || (Case |l|), the three vectors z~", n = 0, 1, 2, in 
(T) are cyclic for the three disjoint representations the original representation 
decomposes into, respectively. For the second case in Table | (Case |), the pair 
z*^, has the same property, as well as the pair z~^, or any pair of the form 
Az*^ + /iz~^, z~^. So far we do not know a single example (for general N, g) where 

B (JC) is not abelian, or is not an algebra. 



6. Reduction of representations as a reduction of N 

Consider an arbitrary element A e 'P{T,\J{N)). Hence both and the genus 
g are arbitrary, but given. We saw that there is an associated wavelet filter ■m^'^\ 
as well as a representation T^"*-' of On acting on L^(T). (As before we work in 
the standard Fourier basis for L'^(T), i.e., e„(z) — z", n E Z.) We identified the 
subspace /C spanned by{e_r,... ,e_i,eo} where r was picked so that reducibility 
of T^"*^ is decided by the fixed-point set of an associated completely positive map 
o-jc given by 

N-l 

(6.1) 4^^(-) = E^/^^(-)^/-'^* 

i=0 

with 

(6.2) V^^^'> ^Vk.T^^\ 

(A) 

While we identified the complete spectral picture of crjj- in the special case g — 2, 
the higher-genus case is not yet entirely understood. But we note that the identity 
from the g = 2 case, 

(6.3) 4^)(So,o) = Ao(yl)£;o,o, 



which is clear from ( 5.63 ), remains true for arbitrary g. Hence our use of cr^'' yields 
reducibility of the representation T^"*^ whenever \o{A) = 1. We recall that if 

(6.4) A{z) = + zyl(i) + • • • + z^-^A'-^-^^ , 

then Xa{A) is defined as the (0, 0)-matrix entry in A'^"'>*A'^°\ 

Remark 6.1. To prove (|6.3|) for the general case, recall that 



N-l N-lg-1 

(6.5) li^^eo = E A, (z^) e, = 5: E A^^e,^,^. 

j=0 j=0 k=0 

Using (|6^), we therefore get 

AT-lg-l 

(6.6) Vl^^eo = P^T^^^^o = E E <^ P^^^^>^n = A^^^eo, 

j=0 k=0 
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and by ( ^l6| ), 

(N-l 
■i=0 

\ 1 0,0 




|eo) (eol 



which is (6.3) 



But it is worth stressing that the argument (6^), f^or the vector eo in /C, does 
not carry over to the other basis vectors e_i, e_2, • • • • This is clear aheady in the 
case g — 2 from ( ^.63 ). More generaUy, for e_i, for example, 



(6.8) 

and therefore 



(6.9) Cr^'^-' (i?_i,_i) — Aq"^£'o,0 + -^TV-l, AT -1-^-1-1 + -^^-2,^-2-^-2, -2 + 

+ '^AT-l.O-'^O,-! + AQjv_iii-1^0 + Ajv_2,Af-l-^-l,-2 + Ajy_;^_jy_2^-2,-1 H , 

where A.^^^ := (A^'^) ., i.e., the («, j) entry in the iV x TV matrix A'^^'^*A^'-\ 



,(00) 



,(00) 



(00) 



Even if A 

(|6j|) will necessarily vanish; see the details below. While 



= 1, that does not imply that all the other coefficients A^-^'' from 



3.10) A^°-'i 



A 



(00) 



A 



the other coefficients in (6.S) such as 



(00) 
(00) 







for aUi,je {0, ...,7V-2}, 



or A, 



(11) 



'^jv-2 JV-2 •'^0 would typically be nonzero 
even if A^''^^^ — 1. So even then, will not be fixed by and there 



is not a natural analogue to (6.3). While if g = 2, then the difference A 

(00) 



(00) 

Af-l,JV-l 



•^^-2 ^-2 ill ^^"^ spectrum of see Table this will generally not be true if 
g > 2. Then the spectral picture for cr^-' is not fully understood. 



However, the following result shows that a general wavelet representation T''^) 

K 



may be reducible because a special point Xo{A) in the spectrum of ct^-* may be 



one. This is a special reduction, however, as we also showed that whenever any 
point A in spectrum ^cr^"*^ satisfies A = 1, we get a possibly different reduction of 
the representation (T(^), ^^(T)). 

Theorem 6.2. Let A e V{T,\]{N)). Then the following jive conditions are 
equivalent. 

(i) MA) = 1- 

(ii) There is a B ^ (^«j)»^j=\ e 7'(T, U(A^ - 1)) such that 

\ 







/ 1 


• 












Sia(z) ••• • 


i?l,Ar-l(-z) 


(6.11) 


A{l)-'A{z) = 













V 


Bn-i.i{z) ■■■ ■ 


Bn-i,n-i{z) 
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I.e., 



(6.12) 



A{l)-'A{z) 



V 



B{z) 



J 



{l)®B{z). 



(iii) After modifying with a \][N)- automorphism of On, the corresponding 
wavelet filter m*-"^' satisfies 



(6.13) 
(iv) 

(v) 



TOq'^''(z) = 1, for all z G T, 



Tj^"*'' Co G Ceo, for all i — 0, . 



4^'*(So,o)eC£;o,o. 



,7V- 1. 



Remark 6.3. Formula (6.13) explains the assertion made in the Introduction 
about reduction in the size of the number N which serves as the scaling of the 
wavelet in question. Once B is identified then there is an {N — l)-wavelet filter 

Ar-2 

(6.14) mPl(z)^ ^S,,+l(z^-l)z^ 
i = l,... ,N -I. 

Proof. (§) ^ (P) : Using the factorization (5.22) in Corollary 5.2, we note that 

(6.15) = (1 - Qg-l) (1 - ■ ■ ■ 



{t-Qi) (a-Qo) (1-Qi)' 



il-Qg-2)il-Qg-l) 



and the V (e V{N)) from ( p. 231 ) is F = A(l), i.e., evaluation of A(z) at z = 1. Let 
£o be the first canonical basis vector in C^. Then (|) states that 

(6.16) ||(l-Qo)(l-Qi)---(l-Q3-i)£o|| - 1. 
Hence ((1 - Qg-i) £o I (1 - Qg-i) to) = 1, where 

(6.17) H:={1- Qg-2) ■■■{!- Qi) | (1 - Qo) | (1 - Qi) • • • (1 - Qg-2) ■ 

Using Schwarz's inequality and induction, we conclude that (Iat — Qj) Eq = Eq, and 
therefore QjEo — for all j = 0, . . . , g — 1. Hence (Iat — Qj + zQj) eq = eo, and 
from ( |5.23D , 

(6.18) A(l)-M(z)eo =£0- 
The same argument yields 

(6.19) {Eo\A{l)-'A{z)ej) =So,j, 
and so (0) follows. 
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Substituting ( |6Al| ) into ( ^.24D , we get 



(6.21) 
and 

(6.22) {V-'A)^Jz) = B,^j{z) ifi,j>l, 

and (|6.13 ) i n (pli[ ) follows. 



(iii) (hj): Using 
(6.23) 7;(^) *e(z) = ^ E ^^fV)C(^"), e e i2(T), 



and introducing the matrix product from ([i^ , we get 



(6.24) 



Using again (2.24), wc thus get 
(6.25) 



* 77 — 



where Vi j denotes the matrix entries of V (:= ^(1)) in \J{N). 
& ^ (B) ^ d) : By ( p3|) , 



(6.26) 

and so if (|^) is assumed, then 



(6.27) 
Since 
(6.28) 

in general, we get 
(6.29) 



^fg) =0 for A: > 0. 



J2Ai,o{z)Ai^Q{z) = 1 

i 

Ao(^) = E^<o^ = l- 



From (3.26), we get 
(6.30) 

and by ( p7| ), 

(6.31) 4^)*(i?o,o) = EEC^-^-' 



fc=0 



fc=0 z=o 



Note the coefficient of i?o,o in (S.31) is Aq (A) — xI^q' ■ For the other coefficients, 
we have 



(6.32) 



(ki) 



.{kk).{ii) 

— ^0,0 ^0,0 
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and 
(6.33) 



^0,0 — K^i^O 



i=0 



Hence, the vanishing of the non-(0, 0) coefficients in ( 3.31) is equivalent to condition 
(3.27), which we already showed is equivalent to both (|y|) and (|). □ 

For ^ e /C, let [^] denote the cyclic subspace in L'^ (T) generated by ^ and 



the operators Tj;^\ and let E';^' [£] be the corresponding projection. 



Corollary 6.4. Suppose condition (|) of Theorem | 
operator W in the commutant of T^'^^ which satisfies 



holds. Then the {unique) 



(6.34) 



isW ^ E^^^ [ea] 



Proof. Let v £ be the vector with coordinates Vi :— A^^. Then (|) holds 

if and only if ||ti|| = 1. In that case, the restriction of r(^) to Ti^'' [eo] is the 
subrepresentation which is induced by the Cuntz state |Cun77| ujy. Specifically, if 
/ = {ii, . . . ,ik), J = {ji, ■ ■ ■ are multi-indices, then it follows from (|^) that 



(6.35) 



J/ J, 7 eo 



■^^k^3l ■■■Vj, {sjS*j). 



Hence [eo] is the closed subspace spanned by eo and |T|'^''eo|. Since 
(6.36) 



N-l 

^A,,,(z^) z^ eP(T,C) 

3=0 



Ti^'' [eo] is contained in the Hardy space 7i+ := span{e„ ; n > 0}. Since /C — 
spanje^r, • ■ ■ , e_i, eo}, it follows that E^^ [eg] satisfies ( |6.34 ). Since fi^"^' (-£-0,0) = 
i?o,Oj the operator W in the commutant of T*^"*^ satisfying (3.34) is unique by 
Theorem Ol and so W = E^_^^ [eo]. □ 



Remark 6.5. The Hilbert space prior to Corollary xA is of the form TY^-* [^] 
for some S, (z IC where 



T 



[Vi) £ 



(6.37) 

for some v 

is unique from 1; up to unitary equivalence, by [ |Cun77 |. It also follows from [ Jor99 
Theorem 6.3] that every £^ £ JC which satisfies (6.37) for some v £ must be a 



Since tOy is a Cuntz state, the corresponding representation 



monomial, i.e., ^ (z) 



for some k £ {0,1, ... , r}. In our special setting, this 



can be proved as follows: If T*^ = ViS, for a unit vector S, £ )C, where J2i l^^l = 
then i = J2^ T,T*^ = Y.^ Or, spelled out, ^ (z) = Y.^ v^^n^ {z) C (z^) . Hence, 

putting TO (z) = 'Y^^Vimi{z), we have ^ (z) — m{z)S^{z^^. Now we may use 
the argument from |Jor99| , page 104] or |BrJo97b, Theorem 3.1] to conclude 
that |to(z)| — 1 and [^ (z)| = 1 for almost all z e T. But both to (z) and ^ (z) 
are polynomials, and we deduce from Lemma 3.1 that they are monomials. If 
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^{z) — z ''it follows that to (z) = ^ (z) ^ (z^) — z^^ Thus a very restrictive 
necessary and sufficient condition for the Cuntz state LUy to occur is that 

(6.38) VwjTO, (z) = z(^-i)'= for some fee {0,1,..., r}. 



The other details of the decompo sitio ns may be spelled out as follows: Let 
H = (T). Note first that Corollary |6.4| holds whenever Eq q is replaced by any 
one-dimensional projection in the fixed-point set of cr in the generality of Theorem 



3.1, since £^o,o (^) — Eq^ (Ti) is then a one-dimensional T^^^^ *-invariant space. If 



e is a general (not necessarily onc-dimcnsional) projection in the fixed-point set. 



it follows from Lemma 3.3 in BJKWOO | that e commutes with all the operators 
prp(A)*p prp{A) p hence the projection E onto [Oi^elC] in H. is in the 
commutant of the representation, and PEP = e. This justifies claims in Remark 
5.11. In the case (7 = 2, which is completely enumerated in | BEJOO| (for = 2) 



and in Table g (for TV > 2), the ranges of the projections in the fixed-point set 
are spanned by subsets of the orthonormal set {eo,e_i, . . . ,e_r}, and hence each 
Hilbert space in the decomposition has the form E'^f^ [^] for a suitable ^ in this set. 
The vector ^ defines a Cuntz state if and only if the corresponding projection in 
B {icy is one-dimensional. Note also that the result por99| . Theorem 6.3] spelled 
out above implies that all one-dimensional projections in B [K-Y are diagonal in the 
standard basis. This gives a partial explanation of these results. 

Other examples of decompositions of (T). Generically, the wavelet 
examples do not have minimal projections e in B {K-Y of dimension one, and even 
if there is one such projection in B {ICY , there may be others that do not have 



dimension one. For illustration, let us repeat some of the examples in |BEJOO|, 



rephrased in the setting of the present paper. Take, for example, N ^ 2 ^ g. 



Then, by (6.38), only the two cases [ ^ and M have the Cuntz-state 



vectors. For A (z) ~ j ' even though Aq (A) = 1, this example has a minimal 

2-dimensional projection e in B{1CY , namely e = -I- E-2,-2- It is minimal 



in the sense that the restriction of T^^^ to E^ — [C'2e/C] is irreducible; see [BEJOO 



(4.52)]. The matrix A{z) ^ ^ ) has Aq {A) = 0, and its representation T^^^ 



^z 0^ 

decomposes into a sum of two irreducibles associated with respective 2-dimensional 
projections e and / in B{K,Y: e = i?o,o + and / = S-2,-2 + -B_3,_3; see 

|BEJOO| , (4.46)]. Hence in this case, H = (T) = [Oze/C] [O2//C] with T^-*) 



restricting to an irreducible representation on each of the two subspaces. A direct 
application of ( ^.15 )-( 5.16| ) in Theorem 5.1 shows that these two subrepresentations 



are disjoint, i.e., inequivalent. While these examples have .9 = 2, the cases g > 2 
entail a richer decomposition structure. 



Note that the result. Theorem 5.1, which is used in analyzing the representa- 



tions T''^^ applies also when g > 2. Moreover, the factorization result. Corollary 



5.2, yields in principle a way of understanding the general case, i.e., arbitrary g and 
N, but unpublished experimentation with examples for g = 3 (i.e., two projections 
in C^) has not so far yielded decomposition structures more general than the above 
mentioned ones. 
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Remark 6.6. The proof of (6.3) shows more generaUy that 

/N-l \ 



(6.39) 



(-Eo.o) 



(0) r(0) 



1=0 



Referring to Theorem |5.1| , wc note then that the necessary and sufficient condition 
for _Eo,o to induce an operator in (T) which intertwines the two representations 
T^-^) and T(^) is (^(") ^ = 1. Moreover, if this holds, then both T^^' and 

T^^^ must satisfy the equivalent conditions in Theorem 3.2, and we must then 
further have v4,-°q — sf^ for all i. Then the intertwining operator W which is 
induced by -Bq.o via a^-^'^^^ (E^.o) — ^'O.o is the one which results from the uniqueness 
of the ^^-representation where w,; = A^^^ = B^j^ . 
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